#ifndef _convective_h_
#define _convective_h_

double calc_r(double i); 

double calc_vr_center_average(sv_t *y, int i, int j, int k);

double calc_vt_center_average(sv_t *y, int i, int j, int k); 

double calc_vz_center_average(sv_t *y, int i, int j, int k); 

double calc_vr_edge_t_average(sv_t *y, int i, int j, int k); 

double calc_vz_edge_t_average(sv_t *y, int i, int j, int k); 

double calc_vr_edge_z_average(sv_t *y, int i, int j, int k); 

double calc_vt_edge_z_average(sv_t *y, int i, int j, int k); 

double calc_vz_edge_r_average(sv_t *y, int i, int j, int k); 

double calc_vt_edge_r_average(sv_t *y, int i, int j, int k); 

double convective_central_r(sv_t *y, int i, int j, int k);

double convective_central_t(sv_t *y, int i, int j, int k);

double convective_central_z(sv_t *y, int i, int j, int k);

#endif

